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We consider a Landau theory of coupled charge and spin-density wave order parameters as a simple 
model for the ordering that has been observed experimentally in the La2Ni04 and La2Cu04 families 
of doped antiferromagnets. The period of the charge-density wave is generically half that of the 
spin-density wave, or equivalently the charges form antiphase domain walls in the antiferromagnetic 
order. A sharp distinction exists between the case in which the ordering is primarily charge driven 
(which produces a sequence of transitions in qualitative agreement with experiment) or spin driven 
(which does not). We also find that stripes with non-coUinear spin order {i.e. spiral phases) are 
possible in a region of the phase diagram where the transition is spin driven; the spiral is circular 
only when there is no charge order, and is otherwise elliptical with an eccentricity proportional to 
the magnitude of the charge order. 



Experiments on the doped lanthanum nickelate pi and 
lanthanum cuprate |y,^ famiUes of materials have estab- 
lished that topological stripe phases are a prominent fea- 
ture of doped antiferromagnets. In the low-temperature 
stripe phase, the doped holes are concentrated in peri- 
odic walls which are simultaneously discommensurations 
in the Neel order. In this paper, the phase diagram of 
coupled charge density wave (CDW) and spin density 
wave (SDW) order parameters will be constructed from 
the Landau theory of phase transitions W. By associ- 
ating the experimental observations with distinct regions 
of the global phase diagram it is possible to gain insight 
into the microscopic mechanism of stripe formation and 
the nature of the spin ordering. It will be shown that 
the experiments are consistent with the suggestion that 
stripes are produced by frustrated phase separation Iq-R] 
and not by a Fermi surface instability g . 

The analysis also addresses the existence of spiral mag- 
netic order found in some theoretical studies of doped an- 
tiferromagnets 0. A circular spiral phase (in which the 
magnitude of the ordered moment is a constant) is only 
possible if there is no accompanying charge order; the 
coupling of the spin and charge order generally produces 
an elliptical spiral phase, with the eccentricity of the el- 
lipse proportional to the magnitude of the charge order 
parameter. In the region of the phase diagram which we 
associate with all known experiments, spiral phases are 
a remote possibility. They may appear only as a third 
transition (yet undetected) at very low temperatures. 

It will be assumed that the Landau free energy func- 
tion depends on only the fundamental Fourier compo- 
nents of the SDW and CDW order parameters. Of 
course, as always, higher harmonics will appear as the 
magnitude of the order increases, but this has no effect 
on the nature of the phase diagram. So as to focus on 



the situation of immediate experimental relevance, we 
will consider the quasi-two dimensional case in which the 
ordering vectors lie in a plane (and, ultimately, in a line), 
but the generalization to other geometries is straightfor- 
ward. (We shall return to the issue of the fluctuation 
effects peculiar to quasi-two-dimensional systems at the 
end.) Also, we consider the case in which the crystal 
has the symmetry of a square lattice, and the ordering 
vectors lie along a symmetry direction, so that there are 
only two inequivalent directions of the ordering vector 
(along the diagonals in the case of the lanthanum nicke- 
lates, and along the vertical and horizontal directions for 
the lanthanum cuprates). This will allow us to extract 
the essential physics in a relatively simple way. The gen- 
eralizations to the various space groups of the structures 
of any given material (which in principle should be used) 
and to allow different ordering vectors is straightforward. 
With these restrictions, the stripe order can be de- 
scribed by the two complex scalars, pj: and pp and the 
complex spin vectors S^- and Sqv , corresponding to charge 

order and spin order respectively. Here, the vectors (fc, 
q) are related to (fc', (f) by rotation through 7r/2, and q 
is measured relative to the magnetic ordering vector Q 
of the undoped system. (It is assumed that Q is unique, 
which requires that Q = —Q, i.e. that 2Q must be equal 
to a reciprocal lattice vector.) 

The most general free energy up to fourth order is con- 
structed by including all terms allowed by symmetry, i.e. 
translation, time reversal, reflection, and spin rotation 
invariance, and the crystal point-group symmetries: 



^ = lfp\p\ 
1 



UM' (1) 

+ [/,|S|4 + L/:,(SxS*)-(S xS*) 



+Ai[(S-S)p* + c.c] +2A2|Sp|p|2. 

Note that it is not necessary to include separate contri- 
butions from wave vectors —q and —k since the charge 
and spin densities are real, and hence p_j^ — pt and 
S_^ = St. Also, there is no separate contribution of the 
form |S • Sp because it may be written as a linear combi- 
nation of the other quartic terms. J-'mt, which is entirely 
quartic, has cross terms coupling the order at {q, k) and 
{(f ,k') ( e.g. yi|Sp||S'p). For simplicity, and because 
it is the case of experimental interest, it will be assumed 
that the interactions in J-int are uniformly repulsive, so 
that unidirectional order is favored (i.e. spin order will 
occur at (for (f, but not both). For example, we do not 
allow Tint to contain terms that favor checkboard order 
p| as an alternative to stripe order. 

The third order term (Ai) coupling spin and charge is 
allowed if and only if 



k^2q. 



(2) 



As discussed previously Cn], this relation is the generic 
reason for the "topological" character of the dopant- 
induced structure, as it implies that the period of the 
spin modulation is twice the period of the charge mod- 
ulation or, in other words, a periodic array of hole lines 
induces an array of antiphase spin domains. 

The free energy in Eq. (1) does not contain umklapp 
terms, in which the sum of wave vectors is equal to a 
reciprocal lattice vector G, which become important in 
the neighborhood of a commensurability. Assuming that 
Eq. (2) is satisfied, the possible umklapp contributions 
up to fourth order are S • S (when 2q — G); p^, S • Sp, 
and (S • S)2 (all when 4g = G); p^ (when 6q = G) and p* 
(when 8g = G). Higher order terms give (weaker) com- 
mensurabilities at smaller wave vectors. As usual [O, 
the system will display commensurate regions separated 
by soliton "discommensurations" when the wave vector 
is close to a commensurate value. 

From now on we will drop the subscripts on the order 
parameters, since each order parameter has a single wave 
vector, as specified above. Also the normalization of the 
order parameters will be chosen so that Up = 1 = Ug- 

Nature of the ordered phases: With only two wave 
vectors (iq), the spin order must be either coUinear or 
coplanar, since a full three dimensional spin texture re- 
quires at least three ordering vectors. We consider the 
coUinear and the non-coUinear cases separately. 

When the spin order is coUinear, the axes and origin 
of coordinates may be chosen so that S — \S\ei and 
p = \p\e^^ , where 9 determines the relative phase of the 
charge and spin density waves. In real space this means 
that the charge and spin density are 



p(r) ~ p — 2\p\ cos(2(f • f— 6 

S(f)e*'5-?^2|S'|eiC0s(q-f) 



(3) 



It is easily seen that the free energy is minimized with 
6» = TT for Ai > and 6' = for Ai < 0. Clearly, for 9 = n, 
the charge density is peaked on the spin domain walls 
where the magnitude of the spin order is zero. Since, 
in all microscopic models studied to date, doping tends 
to depress magnetic order, we expect on general grounds 
that Ai > 0. 

Thus, for coUinear spins, the Landau free energy can 
be expressed in terms of the magnitudes of the order 
parameters as 



Funear{\pl\S\)^-rs\Sf + \S\^ 



(4) 



+ lrM' + \p\' 
-2\X,\\S\'\p\+2X2mp\\ 

For a coplanar spiral phase, the origin of coordi- 
nates and axis of quantization may be chosen so that 
S = |5|[cos(a)ei ±isin(a)e2], which corresponds to a 
real-space spin density 

S{f)e"^-^= (5) 

2|S'|[cos(a) cos(q • r)ei ± sin(a) sin((f • 7^62]. 

Clearly, tan a determines the eccentricity of the ellip- 
tical spiral, and a — 7r/4 corresponds to an ideal spi- 
ral, in which the magnitude of the magnetic order is 
a constant. The minimization of the free energy (with 
the assumption that \S\ ^ 0) with respect to a (in the 
range < a < 7r/4) and 9 can be carried out straight- 
forwardly. If lAillpl/C/^IS-p > 1, the result is a = 0, 
or in other words the coUinear state is recovered. For 
|Ai||p|/C/a;|S'P < 1, the free energy is minimized for 



|^|'cos(2a) = |5ip-|52|'-|Ai||p|/C/, 



(6) 



and as a function of \p\ and \S\ the free energy of the 
spiral state is 



F< 



8piral{\p\,\S\) = ^''^1'^ 



2 + (l-C/,)|5|4 
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(7) 



The spiral phase is limited by the constraint | cos(2q:)| < 
1, which for the simple case of A2 = is satisfied only for 



U, 



2|Ai|(l-f/,) 



2A| 



>0 



(8) 
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FIG. 1. 

The phase diagram shown in Fig. 1 was derived by 
minimizing the free energy in Eq. (1) with respect to p 
and S where, to be concrete, we have shown the locations 
of the phase boundaries for the case A2 = and < Ux < 
1. It exhibits five distinct phases: I) a disordered phase, 
with \p\ = and IS*] = 0; II) a charge ordered phase with 
IpI 9^ and \S\ = 0; III) a collinear stripe ordered phase 
with IpI 7^ and \S\ / 0; IV) an elhptical-spiral phase, 
with IpI f^ and jS'l ^ and < a < n/A as determined 
from Eq. (ph. V) a circular spiral phase with spin order 
but no accompanying charge order. For Ux < 0, the two 
spiral phases are eliminated from the phase diagram, but 
the remaining phase boundaries are unchanged. 

Let us briefly sketch the analysis that leads to this 
phase diagram: 

The Spin Driven Transition to a Collinear Phase: For 
Tp > 0, minimizing Funear with respect to p gives 



IpI 



2Ai|5|2 



0{\Sf), 



(9) 



and, on eliminating p to obtain an effective free energy 
for the spins. 



f: 



•eff 



1 
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linear 



{\S\)^-r^S\^ + {l-^)\S\'' + 0{\S\'). (10) 



Clearly, for rp > 2Af , there is a second order transition 
as a function of r^ from the disordered phase for r^ > 
to the collinear stripe phase for r^ < 0. This transition is 
spin driven; if we assume that Ts oc (T — Tc), we find the 
usual mean-field exponent, \S\ ^ (Tc — T)^/^, while the 
charge modulation, which is parasitic to the spin order, 
grows more slowly, as \p\ ^ jS'p ^ (Tc — T). The second 
order line along Vg = ends at a tricritical point, denoted 
by Ti in Fig. 1, where Vg = and r^ — 2Af , so that the 
coefficients of both the \S\'^ and |S'|'' terms in F^ff{\S\) 
vanish. 



The Transition Driven by Spin-Charge Coupling: Be- 
low the tricritical point, where Vp < 2Af, the transi- 
tion becomes first order, and moves into the quadrant 
in which both Ts > and rp > 0; here, the transition 
is driven by the coupling between spin and charge. The 
precise shape of the first order line depends, in general, 
on A2 and on the values of higher order terms neglected 
in our truncated form of the Landau free energy in Eq. 
(1); the general topology and structure of the phase dia- 
gram shown in Fig. 1 are, however, unchanged by these 
higher order terms. 

The Charge Driven Transitions: For r^ positive and 
sufficiently large, there is a second order line along the 
rp — axis separating the disordered and charge-ordered 
phases. This line ends at a critical end point, Ei in 
the figure, where the first-order line discussed previously 
crosses the axis. In this range of Ts, there is a second 
transition at negative rp from the charge-ordered phase 
to the collinear stripe phase, at which the spin-charge 
coupling finally causes the spin density to order as well. 
Again, one can analyze this transition by first minimiz- 
ing the free energy with respect to p, and then analyzing 
F'^-^-f as a function of 5'. As a result, there is a second 
tricritical point, T2, at which this transition changes from 
first order (an extention of the previously discussed first 
order line) to second order. In either case, the spin order, 
enhances the charge order. In this region of the phase di- 
agram, higher order terms in the Landau free energy can 
have quantatitive ( but, we believe, not qualitative) ef- 
fects on the results. However, to be concrete, it is useful 
to display analytic results obtained with A2 and all higher 
order terms neglected. Specifically, in this case. 



p = po + [Xi/\rp\]\S\' + 0{\S\% 



(11) 



where pq 



is the value of p in the charge- 



ordered phase, T2 is the point r^ — 2\\ and 



-A^ 



and the second order line is given by the expression 
Ts = 2Ai -y/lrpj. The first order line, which is given by the 
expression Tg = Af-|-|rp|, can be located straightforwardly 
once it is realized that on approaching this line form the 
stripe-ordered side, p = Ai/2 and |5p = (l/4)(2Af -r^). 

The Transition to the Spiral Phases: For Ux < 0, it is 
easy to see that a collinear phase always has lower free 
energy than any competitive spiral phase. For Ux > ^, 
the Landau free energy as written is unbounded below, 
so higher order terms must be included in any analysis. 
However, for Ts < and < C/a, < 1, there is the possi- 
blility of a spiral phase. 

Therefore, to complete the phase diagram for r^ < 
and < C/j; < 1, one must compute the minimal 
value of Fapirai (subject to the constraint |cos(2a)| = 
|Ai||p|/f/a;|S'^| < 1) and compare it with the minimal 
value of Funear ■ For simplicity, we first consider the case 
A2 = 0. It is easy then to see that for rp/\\ — 2/Ux > 0, 
Fsptrai is minimized by p = and that Fspirai < Funear ; 
this is region V of the phase diagram, the circular spiral- 
spin phase. At the point where rp/X\~2/Ux changes sign. 



it is clear from Eqs. (6) and (7) that there is a second 
order transition (from region-V to region-IV ) marked 
by the onset of both charge order, and an eUiptical ec- 
centricity to the spin spiral. As Tp decreases, the spiral 
eccentricity gradually increases until it achieves its max- 
imum value (linear polarization). The line cos(2a) = 1 
thus determines the phase boundary between regions IV 
and III of the phase diagram; for A2 = 0, this line is 
simply the parabola determined by Eq. (8). The main 
effect of A2 > (in this region of the phase diagram) is to 
expand the region of the circular spiral phase at the ex- 
pense of the elliptical spiral phase. For A2 > a/(1 — Ux) 
the elliptical spiral phase is completely eliminated. Con- 
versely, of course, A2 < tends to stabilize the elliptical 
spiral phase. 

Further Theoretical Considerations: 
Effect of Higher Harmonics: So far we have considered 
only ordering at the fundamental wave vector, although 
of course higher harmonics are induced below Tj,. It is 
important to verify that these harmonics do not affect 
the stability of the various phases. Slightly below Tc, in 
a spin-ordered phase, the effective free energy to order 
I'S'ngP for the n*'' harmonic is of the form 



Tn 



^n \^nq\ 



[A„ • Snq + C.c] 



(12) 



where r„ > and A„ is a function of Sg' and pj:, if 
there is charge order as well. For a coUinear phase, by 
rotational invariance, A„ oc S, so the induced higher 
harmonics are always parallel to the fundamental. Simi- 
larly, for a spiral phase with unbroken time-reversal sym- 
metry, it is straightforward to see that A„ must lie in 
the ordering plane, so that the planar character of this 
phase is unaffected by higher harmonics. However, if 
time reversal symmetry is broken, then a contribution 
to Aq oc Sq- X S_g- is allowed and the planar phase is 
unstable to the formation of a three-dimensional spiral. 

Goldstone Modes and Fluctuation Effects in CoUinear 
Phases: Landau theory is, of course, mean-field theory, 
so it is important to address the effects of fluctuations 
about the mean-field state. In the absence of commen- 
surability effects, there is a Goldstone mode which re- 
flects the broken translational symmetry associated with 
finite q ordering; commensurability effects, if relevant, 
will produce a gap in this mode, which will be smaller the 
higher the order of the commensurability. In any of the 
coUinear spin-ordered phases, there are two Goldstone 
modes which reflect the broken spin-rotational symme- 
try. Any uniaxial (Ising) anisotropy would produce a 
gap in these modes. 

Thermal fluctuations of these low lying modes may 
not dramatically alter the phase diagram in three di- 
mensions, but in quasi-two dimensional systems, such 
as Lai.6-a;Ndo.4Sra;Cu04, they always destroy the long- 
range order, unless terms that break the symmetry {e.g. 
Ising anisotropy and umklapp scattering ) or inter-plane 
couplings {i.e. three dimensional effects) are included 
in the analysis. However, in many cases, this observa- 
tion is accademic. For instance, the correlation length 



of the two-dimensional spin 1/2 Heisenberg model is 
roughly proportional to exp(J/T) at low temperatures 
|13[ (where J is the exchange integral and a is the lat- 
tice constant). Thus the magnetic correlation length can 
exceed the size of the sample at temperatures of interest 
and, for all practical purposes, the state of the system 
is indistinguishable from long-range order. The effects 
of disorder are also potentially dramatic in two dimen- 
sions where the density wave order will generally break 
up into Lee- Rice domains p4[ , but, again, if the disor- 
der is sufficiently weak, this may be of largely accademic 
interest. 

Goldstone Modes of the Spiral Phase: The elliptical 
spiral phase has, in addition to the one Goldstone mode 
associated with broken translational symmetry, and the 
two Goldstone modes associated with rotations of the 
principal axis of the elipse, an additional Goldstone mode 
associated with rotations about the principal axis. Ising 
anisotropy will open a gap in two of the spin-related 
Goldstone modes, but will leave the third one gapless; 
it requires XYZ anisotropy to fully gap the spin-related 
Goldstone modes. Thus, the presence of a third Gold- 
stone mode, in particular one with an anomalously small 
gap, can be used as a diagnostic for the presence of ellip- 
tical spin order. 

In the circular spiral phase, rotations in the plane of 
the spiral arc equivalent to a translation, so the number 
of Goldstone modes is the same as in the coUinear phase. 
However, the mode related to translational invariance has 
a spin component, and thus is relatively insensitive to 
commensurability effects and to the effects of disorder. 

Relation to experimental results: Several fea- 
tures of the experiments on the cuprates and nickelates 
may be discussed in terms of the Landau theory analy- 
sis. Both Lai.6-xNdo.4Sra;Cu04 B and doped nickelates 
|l| show an onset of charge order prior to spin order. 
Lai.6-a;Ndo.4Sra;Cu04 undergoes two continuous transi- 
tions (onset of charge order and later spin order), and 
therefore, can be associated with a path on the phase 
diagram that goes from phase-I to phase-II, and then 
to phase-Ill through second order transitions. Several 
La2Ni04+5 samples |l| clearly show a first order transi- 
tion associated with the onset of spin order in the tran- 
sition from phase-II (|p| ^ , \S\ = 0) to the full stripe 
order, phase-Ill. Thus, these experiments can be associ- 
ated with a path on the phase diagram that goes through 
the first order transition line between point Ei and the 
tricritical point T2. In both cases, it is clear that the tran- 
sitions are charge driven rather than spin driven. Note 
that in this region of the phase diagram there is very lit- 
tle prospect of having a third phase transition to a spiral 
state, unless it occurs at much lower temperatures. 

The above discussion has said nothing about the par- 
ticular value of the wave vector k = 2q. In the Lan- 
dau theory, the value of q can only depend on the (non- 
specific) g-dependence of the coefficients of the various 
terms in the free energy, or on commensurability. Ex- 
perimentally, in the cuprates ||,0l, the ordered stripes 



in Lai.6-2;Ndo.4Sr2:Cu04 and the fluctuating stripes in 
La2_a;Sr2;Cu04 have q w x for x < 1/8, which corre- 
sponds to one hole per two Cu along the stripe. For 
X > 1/8, q « 1/8 in La2_xSrj:Cu04 and somewhat larger 
in Lai.6-2;Ndo.4Sra;Cu04. (Here, q is measured in units 
of 27r/a.) This suggests that the value of q is influenced 
by the internal structure of the stripe for x < 1/8, and 
by the p^ umklapp term for x > 1/8. Furthermore, 
the temperature dependence of q in the ordered phase of 
Lai.g_a;Ndo.4Sr2;Cu04 is not strong, which suggests that 
the amplitude of the stripe is well established at the or- 
dering temperature, and that the transition is produced 
by phase ordering, as expected for quasi two-dimensional 
systems. In all respects, the value of the stripe wave vec- 
tor q is determined by charge dynamics, rather than spin 
dynamics. 

The conclusion of this analysis that stripes are charge 
driven rather than spin driven supports the idea that 
the driving force is Coulomb-frustrated phase separation 
P,p| driven by the hole dynamics. In order to minimize 
their kinetic energy, the holes attempt to separate into 
hole-rich regions and regions with significant local an- 
tiferromagnetic correlations. However phase separation 
is frustrated by the long-range Coulomb interaction be- 
tween the holes, and the compromise between these two 
forces is the (charge-driven) stripe phase. The antiphase 
ordering of the spin domains is a consequence of trans- 
verse stripe fluctuations; it ensures that, on average, ad- 
jacent spins are antiparallel, whatever the location of the 
stripe. 

The alternative mechanism for stripe formation is a 
Fermi surface instability ||^ (due to nesting) , in which the 
spins form antiphase domains that are stabilized by holes 
bound to the domain walls. In this scenario, spin and 
charge order together or charge stripe order follows after 
spin order, contrary to experiment. These microscopic 
theories also predicted the possibility of spiral phases ||] 
(though not that they must be elliptical). We conclude 
that theories based on a Fermi surface instability may be 
relevant in the spin-driven region of the phase diagram 
(rp > 0) but not for any of the materials of interest, 
which all show a charge driven sequence of transitions. 

Indeed there are more general reasons to believe that 
the ordered phases in the high-Tc cuprates should not 
arise from a Fermi surface instability. It has been argued 
|p2| that these materials belong to a class of "bad met- 
als" , in which there are no quasiparicles, and therefore 
no Fermi surface or Fermi surface instabilities. 
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